Walks in the quarter plane: Kreweras' algebraic model 
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Abstract 

We consider planar lattice walks that start from (0,0), remain in the first quadrant i,j > 0, 
and are made of three types of steps: North-East, West and South. These walks are known to 
have remarkable cnumcrative and probabilistic properties: 

- they are counted by nice numbers (Kreweras 1965), 

- the generating function of these numbers is algebraic (Gessel 1986), 

- the stationary distribution of the corresponding Markov chain in the quadrant has an 
algebraic probability generating function (Flatto and Halm 1984). 

These results are not well understood, and have been established via complicated proofs. Here 
we give a uniform derivation of all of them, which is more elementary that those previously 
published. We then go further by computing the full law of the Markov chain. This helps to 
delimit the border of algebraicity: the associated probability generating function is no longer 
algebraic, unless a diagonal symmetry holds. 

Our proofs are based on the solution of certain functional equations, which are very simple 
to establish. Finding purely combinatorial proofs remains an open problem. 

1 Introduction 

Let us begin with a very simple combinatorial statement: the number of planar lattice walks that 
start and end at (0,0), consist of 3n steps that can be North-East, South, or West, and always 
remain in the nonnegative quadrant i, j > is 



a(3n) 



3n 



(n + l)(2n + 1) V n 

An example of such a walk is given in Figure ^ This result, first proved by Kreweras in 1965 |24| . 
is rather intriguing, for at least two reasons. 

First, this simple looking formula has no simple proof. If we consider instead the more traditional 
square lattice walks (consisting of North, South, East and West steps), then there exists a nice 
formula too: the number of 2n-step walks, starting and ending at the origin and confined in the 
first quadrant, is 

1 /2n + 2 N " 

[ U ' ~ (2n + l)(2n + 4) \ n + l 

* Partially supported by the European Community IHRP Program, within the Research Training Network "Alge- 
braic Combinatorics in Europe", grant HPRN-CT-2001-00272. 
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Figure 1: Kreweras' walks in a quadrant. 

But the latter formula can be proved in a few lines (count first the number of such walks having 
2m horizontal steps, and then sum over all values of m) and admits even a direct combinatorial 
explanation [111 \21\ . No similar derivation exists for the numbers a(3n). 

The second fact that makes the numbers a(3n) intriguing is that their generating function, that 
is, the power series A(t) = ^2 n a(3n)t n , is algebraic. This means that it satisfies a polynomial 
equation P(t,A(t)) = where P is a non-trivial bivariate polynomial with rational coefficients. 
For combinatorialists, objects that have an algebraic generating function are really special: this 
property suggests that one should be able to factor them into smaller objects of the same type, and 
then translate this factorization into a polynomial equation (or a system of polynomial equations) 
defining the generating function. Let us take an example: it is known that for any (finite) set of 
steps, the walks confined in the upper half-plane have an algebraic generating function. There is a 
clear combinatorial understanding of this property: the key idea is to factor the walk at the first 
time it returns to the x-axis. It is still an open problem to find an explanation of this type for the 
algebraicity of the series A(t). Let us underline that not all walks in the quadrant have an algebraic 
generating function: the generating function for the numbers b(2n) is transcendental (see [5] for a 
stronger result). 

A natural question — at least for a computer scientist — is whether the set of words on 
the alphabet {a, b, c} that naturally encode Kreweras' walks forms an algebraic (or: context-free) 
language (223- These words contain as many a's as o's, as many a's as c's, and each of their prefixes 
contains no more 6's than a's, and no more c's than a's. Using the pumping lemma [221 Theorem 
4.7], one can prove that this language is not algebraic. Moreover, the words satisfying only the 
second condition above, which encode walks ending anywhere in the quadrant, do not form an 
algebraic language either j3J. However, we shall see that their generating function is algebraic. 

Then how does one prove Kreweras' formula? In his original paper, Kreweras considered n- 
step walks in the quadrant going from (0,0) to A step by step construction of these walks 

gives an obvious recurrence relation for their number, denoted below Oj,(n). Kreweras solved this 
recursion. His proof involves guessing a substantial part of the solution, and then proving several 
hypergeometric identities. The latter part was then simplified by Niederhausen |29l 130] . A different 
proof, due to Gessel, also requires guessing the bivariate generating function of walks ending on 
the x-axis, and then verifying that it satisfies a certain functional equation |20| . 

On the probabilistic side, in the early 70's Malyshev began to address the very general problem of 
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computing the stationary distribution of discrete homogeneous Markov chains in the quadrant [27] . 
Several instances of this question actually correspond to finding the equilibrium behaviour of double- 
queue processes EH EEH E2 This work culminated in 1999 with a book that is entirely devoted 
to solving this problem in the case of unit increments (Fayolle, Iasnogorodski and Malyshev |14j). 
The techniques used in this book are far from elementary, involving sophisticated complex analysis, 
Riemann surfaces and boundary value problems. Solutions are often expressed in terms of elliptic 
functions. The book lists a number of cases in which the stationary distribution has a rational 
generating function, and mentions exactly one case (actually due to Flatto and Hahn JH]) where 
this generating function is algebraic. Not surprisingly, the set of increments of this random walk is 
the same as in Kreweras' problem 1 . 

Hence the following question: what is so special with this set of three steps? Could one find a 
single argument that proves both the algebraicity of the generating function that counts these walks 
and the algebraicity of the generating function for the stationary distribution of the corresponding 
Markov chain? 

This is the question we answer — positively — in this paper. For both the combinatorial 
problem and the probabilistic one, it is very easy to establish a functional equation defining the 
generating function. We solve both equations using the same approach. The only difference is 
that we are dealing with formal power series in the first problem, but with analytic functions 
in the second one. Our solution is constructive (we do not have to guess anything) and more 
elementary than the previously published ones. In particular, we always remain in the (small) 
world of algebraic functions, and do not need to introduce elliptic functions. The key to our 
approach is the combination of the kernel method (which is also central in JH] or ^1]) with a 
special property of the kernel of the equations we consider. Moreover, after having solved the 
counting problem (Section^ and the probabilistic one (Section |3J), we combine both viewpoints 
and compute explicitly the full law of the Markov chain (Section 0J). This actually marks the 
end of algebraicity: the probability generating function is transcendental, unless the transition 
probabilities satisfy a diagonal symmetry. Still, this generating function belongs to the nice class 
of D-finite (or: holonomic) series, which is defined below. 

Obviously, since this paper aims at explaining why a specific set of steps has such special 
properties, it cannot compete in generality with the strength of the machinery developped in |14j . 
Still, it is natural to ask how far our approach could be generalized. It is actually while fighting 
with Kreweras' walks that I discovered it. But it turns out that other applications of this approach 
were published before I was able to complete the present paper. In particular, most of the results 
in Section El are already reported in some conference proceedings together with a general holonomy 
criterion for the enumeration of walks in the quadrant [3]- More recently, the same ideas were 
applied to certain counting problems on permutations [£]. Four new equations were thus solved in 
a uniform, elementary way. Their solutions are usually transcendental, but holonomic, and can be 
expressed as integrals of algebraic (quadratic) functions. I have not tried to attack other stationary 
distribution examples. But it is likely that the method presented here and in [5j can be applied to 
solve explicitly (and in an elementary way) at least certain specific examples. 

Let us conclude this section by giving some definitions and notation on formal power series. 
Given a ring L and k indeterminates xi, . . . , Xf., we denote by L[xi, . . . , Xk] the ring of polynomials 

X A (partial) algebraicity criterion is actually given in [141 Thm. 4.3.1 and 4.3.6], but it is only illustrated by 
Kreweras' example. 
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in x\, . . . , xt with coefficients in L. We denote by L[[sci, . . . , Xk]] the ring of formal power series in 
the Xi, that is, of formal sums 

E a(n l ,...,n k )xl 1 ■■■xl k , (1) 
n\~>0,...,n^>0 

where a(n\, . . . , n^) 6 L. A Laurent polynomial in the Xj is a polynomial in both the Xj and the 
Xi = 1/xj. A Laurent series in the Xi is a series of the form Q in which the summation runs over 
n i > m i for all i, with m, in Z. For F G IL[[£]] , we denote by [t n ]F the coefficient of t n in F(t). If -F 
is a formal series in t whose coefficients are Laurent series in x, we denote by F + the positive part 
of F in x, that is, 

F = E * n E ^ n ) x * F+ = E tn E ^( n ) xi - ( 2 ) 

n>0 ieZ n>0 i>0 

We define similarly the negative, nonnegative and nonpositive parts of F. 

Assume, from now on, that L is a field. We denote by L(xi, . . . , xt) the field of rational functions 
in x\, . . . , Xk with coefficients in L. A series F in L[[xi, . . . , Xk]} is algebraic if there exists a non- 
trivial polynomial P with coefficients in L such that P{F,x\, . . . ,xj.) = 0. The sum and product 
of algebraic series is algebraic. The series F is D-finite if the partial derivatives of F span a finite 
dimensional vector space over the field L(xi, . . . , x^); see (31] for the one- variable case, and |2S1l2f)] 
otherwise. In other words, for 1 < i < k, the series F satisfies a non-trivial partial differential 
equation of the form 

where P^i is a polynomial in the Xj. Any algebraic series is D-finite. The sum and product of 
D-finite series are D-finite. The specializations of a D-finite series (obtained by giving values from 
L to some of the variables) are D-finite, if well-defined. Finally, if F is D-finite, then any diagonal of 
F is also D-finite (the diagonal of F in x\ and X2 is obtained by keeping only those monomials 
for which the exponents of x\ and xi are equal). We shall use the following consequence of the 
proof of this result: if F{t,x) € L[x,x][[i]] is algebraic (with x = 1/x), then the positive part of F 
in x is D-finite, as well as the coefficient of x % in this series, for all i. 

2 Enumeration: the number of walks 

Consider walks that start from (0, 0), consist of South, West and North-East steps, and always stay 
in the first quadrant (Figure Let ajj (n) be the number of n-step walks of this type ending at 
We denote by Q(x,y;t) the complete generating function of these walks: 

Q(x,y;t):= a ilj {n)x i y i t n . 

i,j,n>0 

We can construct these walks recursively, by starting from (0, 0) and adding a step at each time. 
This gives the equation: 

Q(x,y;t) = 1 + t ( - + - +xy) Q(x,y;t) - - Q(x,0;t) - - Q(0,y;t). 
\x y J y x 
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The first term in the right-hand side encodes the empty walk, reduced to the point (0, 0). The next 
term shows the three possible ways one can add a step at the end of a walk. However, one should 
not add a South step to a walk that ends on the x-axis: the third term subtracts the contribution 
of this forbidden move, and the last term takes care of the symmetric case. Equivalently, 

(xy -t{x + y + x 2 y 2 )) Q(x, y; t) = xy - xtQ(x, 0; t) - ytQ(0, y; t). (3) 

We shall often denote Q(x,y;t) by Q(x,y) for short. Let us also denote the series xtQ(x,0;t) by 
R(x;t) or even R(x). Using the symmetry of the problem in x and y, the above equation becomes: 

(xy -t{x + y + x 2 y 2 )) Q(x, y) = xy - R(x) - R{y). (4) 

Equation © is equivalent to a recurrence relation defining the numbers aij(n) inductively with 
respect to n. Hence, it defines completely the series Q(x,y;t). Still, the characterization we have 
in mind is of a different nature: 

Theorem 1 (The number of walks) Let W = W(t) be the power series in t defined by 

W = t(2 + W 3 ). 

Then the generating function of Kreweras' walks ending on the x-axis is 

1/1 1 f I 1 



Q(x,0;t) = — [- 777-- V 7 ! - xW 2 

^ y ' ' tx \2t x \W xj 

Consequently, the length generating function of walks ending at (i, 0) is 

where Ci = (^*) /(i + 1) is the i-th Catalan number. The Lagrange inversion formula gives the 
number of such walks of length 3n + 2i as 

4 n (2i + l) {2i\{?>n + 2i 

a ifi {3n + 2i) - 



(n + i + l)(2n + 2i + 1) \ i J \ n 

The aim of this section is to derive Theorem ^ from the functional equation (j3J) . Note that the 
complete generating function Q(x,y) can the be recovered using 0: 



_ (l/W - x)y/l - xVP- + (l/W -y)y/l- yW 2 1 
xy — t(x + y + x z y z ) xyt 

with x = 1/x. For walks ending on the diagonal, we shall also obtain a nice generating function: 

Theorem 2 (Walks ending on the diagonal) Let W = W(t) be defined as above. Then the 
generating function of Kreweras ' walks ending on the diagonal, defined by 



Qd{x;t) := ^2 a i ^ i (n)x l t n . 



i,n>0 

satisfies 



„ , \ W-x 
tQJx; t) = . = + x. 

V 7 ! - xW{l + Wyl) + x 2 W 2 /4 
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The expression of Qd becomes a bit simpler if we express it in terms of the unique power series 
Z = Z(t) satisfying Z = 1 + 4t 3 Z 3 . Then W = 2tZ and 

_ , « ItZ-x 

tUAx\ t) = ; + X. 

^/l-xtZ{l + Z) + x 2 t 2 Z 2 

The last formula of Theorem ^ is due to Kreweras [2U • He also gave a closed form expression for 
the number of walks containing exactly p West steps, q South steps, and r North-East steps, that 
is, for walks of length n = p + q + r ending at = (r — p,r — q). This expression is a double 

summation, with alternating signs. We have not found anything simpler. 

2.1 The obstinate kernel method 

The kernel method is basically the only tool we have to attack Eq. (jlj). This method has been 
around since, at least, the 70's, and is currently the subject of a certain revival (see HS1 Ex. 2.2.1.4 
and 2.2.1.11] and ^2] f° r early uses of the method, and ^ Ej for more recent combinatorial 
applications). It consists in coupling the variables x and y so as to cancel the kernel K(x,y) = 
xy — t(x + y + x 2 y 2 ) (which is the coefficient of Q{x,y) in @). This should give the "missing" 
information about the series R{x). 

As a polynomial in y, this kernel has two roots 

T , , . l — tX—\/{\— tx) 2 — At 2 X 9 n N 

Y (x) = y\ > = t + xt 2 + 0{t 3 ), 

2tx 



. . 1 — tx + v(T— tx) 2 — 4t 2 x x 9 9 „, n N 

YAx) = — = x 2 - t - xt 2 + 0(t 3 ). 

w 2tx t v 1 

The elementary symmetric functions of the Yi are 

Y + Yi = ~ - x 2 and Y^Yi = x. (5) 

The fact that they are polynomials in x = 1/x will play a very important role below. 

Only the first root can be substituted for y in @ (the term Q(x,Y\;t) is not a well-defined 
power series in t, because of the negative power of t that occurs in Y\). We thus obtain a functional 
equation for R(x): 

R(x) + R(Y ) = xY . (6) 

It is not hard to see that this equation — once restated in terms of Q(x,0) — defines uniquely 
Q(x, 0; t) as a formal power series in t with polynomial coefficients in x. Equation © is the standard 
result of the kernel method. 

Still, we want to apply here the obstinate kernel method. That is, we shall not content ourselves 
with Eq. ©, but we shall go on producing pairs (X, Y) that cancel the kernel and use the infor- 
mation they provide on the series R(x). This obstinacy was inspired by the book ^1] by Fayolle, 
Iasnogorodski and Malyshev, and more precisely by Section 2.4 of this book, where one possible 
way to obtain such pairs is described (even though the analytic context is different). We give here 
an alternative construction. 

Let (X, Y) (0, 0) be a pair of Laurent series in t with coefficients in some field such that 
K(X, Y) = 0. Recall that, as a function of y, the polynomial K(x, y) is quadratic. Thus let Y' be 
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the other solution of the equation K(X,y) = 0. We define the function by ^f(X,Y) = (X,Y'). 
For instance, if (X,Y) is the pair (x,Yq), then ^(X, Y) = (x,Y\). Similarly, we define Y) = 
(X',Y), where X' is the other solution of K(x,Y) = 0. Note that and ^ are involutions and 
that, in view of (JHJ), X' = Y' = (XF)" 1 . In particular, $>(x,Yq) = (Yi,Yq). Let us examine the 
iterated action of <3? and Vl/ on the pair (x,Yq): We obtain the diagram of Figure |2 




(Y ,x) 



Figure 2: The orbit of (x,Yq) under the action of $ and The framed pairs can be substituted 
for (x, y) in the functional equation. 

All these pairs of power series cancel the kernel. We have already seen that the pair (x, Yq,) can 
be substituted for (x,y) in Eq. (@J). This is also true, but less obvious, for the pair (Yq,Yi): indeed, 
if we write 

Q(x,y;t) = ]T t k+l+m x k -"y k - m a k -^- m {k + l + m), 

fc>max(£,m) 

t k+2m (xy) k - e (t y y- m a k - e , k - m (k + i + m), 

fc>max(£,m) 

and note that YqY\ = x, while tY\ = x + 0(t), then we see that Q(Yq, Y±;t) is a well-defined power 
series in t, with coefficients in Q[x,x]. The same argument shows that R(Y\) = tYiQ(0,Yi;t) is 
also well-defined. Thus the two pairs than can be substituted for (x, y) in the functional equation 
give us two equations for the unknown series R(x): 

f R(x)+R(Y ) = xY , 

\ R(Y )+R(Y 1 ) = Y Y 1= x. {l> 

Remark. Let p, q, r be three nonnegative numbers such that p + q + r = 1. Take x = (pr) 1 / 3 ^" 2 / 3 , 
y = (gr) 1//3 p _2//3 , and t = (pqr) 1 / 3 . Then K(x,y;t) = 0, so that R(x) + R(y) = xy. This equation 
can be given a probabilistic interpretation by considering random walks that make a North-East 
step with (small) probability r and a West (resp. South) step with probability p (resp. q). This 
probabilistic argument, and the equation it implies, is the starting point in Gessel's solution of 
Kreweras problem j^UJ Eq. (21)]. 

2.2 Symmetric functions of Yq and Y\ 

After the kernel method, the next tool in our approach is the extraction of the positive part of 
power series, defined by J5J). This is where the values of the symmetric functions of Yq and Y\ 
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become crucial: the fact that they only involve negative powers of x (see ©) will simplify the 
extraction of the positive part of certain equations. 

Lemma 3 Let F(u, v; t) be a Laurent series in t with coefficients in C[u, v], symmetric in u and v. 
That is, F(u,v;t) = F(v,u;t). Then the series F(Yo,Yi;t), if well-defined, is a Laurent series in t 
with polynomial coefficients in x. Moreover, the constant term of this series, taken with respect to 
x, is F(0,0;t). 

Proof. By linearity, it suffices to check this when F is simply a symmetric polynomial in u and v. 
But then it is a polynomial in u + v and uv with complex coefficients. The result follows, thanks 
to ©. 

■ 

We now want to form a symmetric function of Yq and Y\ , starting from the equations (J7J) . The 
first one reads 

R(Y )-xY = -R(x). 
By combining both equations, we obtain the companion expression: 

JR(Yi) - xYi = R(x) + 2x- l/t. 

Taking the difference 2 and dividing by Yq — Y% gives 

R(Y )-R(Y 1 ) 2R(x) + 2x-l/t 

77 77 X = tX . , (8) 

Y - Y 1 ^/A(x) 

where A(x) = (1 — tx) 2 — 4t 2 x is the discriminant that occurs in both Yq and Y\. 

As a Laurent polynomial in x, A(x) has three roots. Two of them, say Xq and X±, are formal 
power series in y/t; the other is a Laurent series in t (for generalities on the roots of a polynomial 
over C(t), see jHSl Chapter 6]). The coefficients of these series can be computed inductively: 

1287 

X = t + 2t 2 Vi + 6t* + 21t 5 Vi + 80t 7 + —^-t 8 Vt^ 

Xt = t-2t 2 Vi + 6t 4 -21t 5 Vt + 80t 7 -^^t 8 \/i + --- 

X 2 = — ;r - 2t - 12t 4 - 160t 7 - 2688t 10 - 50688t 13 + • ■ ■ 
At 2 



Hence A(x) factors as 
with 



A(x) = A A + (x)A_(x) 



A = 4t 2 X 2 , A + (x) = l-x/X 2 , A-(x) = (l-xX )(l-xX 1 ). 

Note that Ao, A+(x) and A_(x) are power series in t with constant term 1. Moreover, Ao has its 
coefficients in Q, while A + (x) has its coefficients in Q[x], and A_(x) has its coefficients in Q[x], 
This is an instance of the "canonical factorization" of power series of Q[x, x] [[t]], which has already 

2 An alternative derivation of Kreweras' result, obtained by considering the product (R(Yo) — xYo)(R(Yi) — xYi), 
is presented in 
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proved useful in several path enumeration problems |19l |U |U] . Going back to (JSJ) , and multiplying 
through by \J A_(x), one obtains 



( R(Y ) - R(Y{) \ _ 2xR(x) + 2-x/t 
V Y Q -Y 1 X ) ~ ^A^T~ ' 



Both sides of this identity are power series in t with coefficients in Q[sc, sc]. But the right-hand side 
only contains nonnegative powers of x, while the left-hand side, except for a term —x, only contains 
nonpositive powers of x (in view of Lemma EJ). Extracting the positive part of the above equation 
thus gives 

t ( 2xR(x) + 2-x/t \ 



-x 



^0 V V^H 



The expression of Q(x,0) announced in Theorem ^ follows, given that X2 = 1/W 2 and R(x) = 
xtQ(x,0). The expansion of Q(x,0) in x is straightforward, using 1 — y/1 — At = 2t J2 n >o Cnt n - 
The value of ai t o(3n + 2i) follows using the Lagrange inversion formula [331 p. 38]. 



2.3 The algebraic kernel method 

We present in this section another proof of Theorem ^ based on a variation of the kernel method. 
This variation does not require to cancel the kernel, but, instead, builds on one of its algebraic 
properties. This variant has some drawbacks — since the kernel is not zero, we are handling bigger 
equations — but it also has some advantages. In particular, we obtain at some point an equation 
that is the counterpart of (JSJ), but in which it is obvious that the left-hand side is nonpositive in 
x. This will be helpful in the next section, where we handle analytic functions rather than power 
series. Finally, this variant of the kernel method provides a proof of Theorem [2 
Let us return to the original equation or, equivalently, to 

xyK r (x, y)Q{x, y) = xy - R(x) - R(y), 

where K r (x,y) = 1 — t(x + y + xy) is the rational version of the kernel K. The fact that the diagram 
of Figure |21 is nice actually stems from an invariance property of K r : 

K r (x,y) = K r (xy,y) = K r (x,xy) = K r . 

Applying iteratively the (involutive) transformations $ : (x, y) 1— ► (xy, y) and ^ : (x, y) 1— ► (x, xy) 
gives the following set of pairs, on which K r takes the same value: 

Note that the diagram of Figure 12 is the specialization of the above one to the case y = Yq. 

Now, all pairs of the above diagram can be substituted for (x, y) in the functional equation: 
the resulting series are power series in t with coefficients in Q[x,x, y, y\. This gives no less than 
three equations: 



xy 


K r 


Q(x,y) 


= xy 


- R(x) 


- R(y), 


X 


K r 


Q(xy,y) 


= X 


- R{xy) 


R(y), 


y 


K r 


Q{x,xy) 


= y 


- R(x) 


- R{xy) 
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We sum the first and third equations, and subtract the second one, so as to keep R(x) as the only 
unknown function on the right-hand side: 

K r (xyQ(x, y) - xQ(xy, y) + yQ(x, xy)^j = xy - x + y - 2R(x) = - — - — - - 2x — 2R(x). 



Equivalently, 



1 1/1 

xyQ{x, y) - xQ(xy, y) + yQ(x, xy) + - = — ( - - 2x - 2R(x) ) . (9) 



The kernel K(x,y) factors as —tx 2 (y — Yo)(y — Y\). Converting 1/K into partial fractions of y 
yields the following expression for the reciprocal of the (rational) kernel K r : 

,n>0 n>l 



K r y/A(x)\l-yY 1-y/Yy J y/Mx) 



Note that this expansion is valid in the set of formal power series in t with coefficients in Q[sc, x, y, y] . 
Let us extract in © the constant term in y: the series xyQ(x,y) and yQ(x,xy) do not contribute, 
and we obtain 

n . . 1 \lt-2x - 2R(x) 
-xQ d (x) + - = J-—j==-±-L 

where the series Qd is the diagonal of Q(x, y), and counts walks ending on the diagonal. The above 
equation should be compared to (JHJ): basically, both equations are equivalent, but their negative 
parts on the left-hand side are written in two different ways. We now proceed as above, using the 
canonical factorization of A(x), which gives 

1 \ 1/t -2x- 2R(x) 



VAl(E)(j-xQ d (x)^ 



VA A+(x) 



Extracting the nonnegative part gives, as before, the value of R(x), and Theorem^ Extracting 
the negative part gives 

'1 „ , s\ 1 2x 



/AT 
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Recall that A = M 2 X 2 = 4i 2 jW 2 and A_(x) = (1 - xX )(l - xX{), where X and Xi are the 
two "small" roots of A(x). We can express their elementary symmetric functions in terms of the 
third root, X 2 = l/W 2 . This gives 

A_(x) = 1 - xW{l + W 3 /A) + x 2 W 2 /4, 

and this provides the expression of Qd(x) given in Theorem |2j 

■ 

3 Probability: a Markov chain and its stationary distribution 



A 




Figure 3: The transition probabilities. 

We consider a Markov chain on the quadrant, whose transition probabilities T(i, j;k,l) are 
schematized in Figure |21 More precisely, for i > and j > 0, the probability of going from to 
(k,£) is 

{p if k = % — 1 and £ = j 
q if k = i and i = j - 1 

r if k = i + 1 and £ = j + 1 

where p, q, r are three positive real numbers summing to 1. When the point (i, j) lies on the border 
of the quadrant, the transition probabilities are modified as follows: for i > 0, 

T u n- k £) - / P ' ifk = i ~ 1 and £ = 

i\t,u,K,i) - | r , ■ l£k = i + 1 andi= i 

and for j > 0, 

Tf0 7 -. fc ^_ / 9" iffc = and^ = i-l 

MU,J ' '^~\r" iffc = l and£ = j + l 

where p', r', q", r" are positive numbers such that p'+r' = q"+r" = 1. Finally, we take T(0, 0; 1, 1) = 
1. Note that this chain is irreducible (all states communicate) and has period 3. 

A probability distribution (pij)ij>o is stationary for the above transition if for all k,£>0, 

Pk,e = ^2pi,jT(i,j;k,£). 

i,3 
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Our objective is to find the stationary distribution of the above transition, when it exists. It is 
customary to encode a distribution by its probability generating function 

i,3>0 

but is is more convenient here to split H(x,y) into four parts: firstly, po,Oi an d then the three 
following generating functions: 

i,3>l »>1 j>l 

Then the distribution (j*ij)i,j>o is stationary if and only if 

(1 — px — qy — rxy)P(x, y) + (1 — p'x — r'xy)P\{x) + (1 — q"y — r" xy)P2{y) + (1 — xy)po ) o = 0. (10) 

Note that the numbers pij have to sum to 1: hence the above series are absolutely convergent for 
| a; | < 1 and \y\ < 1, and define analytic functions for \x\ < 1, \y\ < 1. Moreover, 

Po,o + A(l)+ J P 2 (l) + i 3 (l,l) = 1. (11) 



3.1 The main results 

The stationary distribution of this Markov chain was computed in ^3] m the case where the 
transition probabilities are related by 

V - = t, = a and q - = ^-:=p. (12) 

Equivalently, 

/ V / r n Q // r 
p = , r = , q = , r = . 

p+r p+r q+r q+r 

It is known that this chain has a stationary distribution if and only if r < min(p, q) (see |281 115j for 
general results on Markov chains in the quadrant). It will be shown in Lemma that the condition 
is necessary, and it will follow from the results of Section^ where we compute the law of the chain, 
that it is sufficient. For the moment, we rely on the general results of |281 115j . 
Under the conditions (|12|) . which we assume to hold in this section, we have 

1 — p'x — r'xy = — (1 — px — qy — rxy + q(y — 1)) , 
r 

r" 

1 — l"y ~ r"xy = — (1 — px — qy — rxy + p(x — 1)) , 
r 

1 — xy = — (1 — px, — qy — rxy + p(x — 1) + q(y — 1)) 
r 

so that the functional equation (|1U|) can be nicely rewritten as 

(1-px-qy- rxy)Q(x, y) = q(l - y)Q(x, 0) + p(l - x)Q(0, y) (13) 

with 

Q(x,y)=p 0fi + r'P 1 (x) + r"P 2 (y) + rP(x,y). (14) 
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This equation was first met by Flatto and Hahn ^S] in their study of a system of two parallel queues 
with two demands (with continuous time). They solved this equation using non-trivial complex 
analysis, multivalued analytic functions, and a parametrization of the kernel by elliptic functions 
— to end up with an algebraic solution Q(x, y). We shall rederive their result in a more elementary 
way, and state it in a more symmetric fashion. 

Theorem 4 (Solution of Flatto and Hahn's equation) Assume r < min(p, q). There exists, 
up to a multiplicative constant, a unique solution of (jl3|) that is analytic in {\x\, \y\ < 1} and whose 
series expansion converges for \x\, \y\ < 1. This solution satisfies: 

^•°) = 71 TVT r r ('fl-^-/i3^5-g/',-J-i./i-, 

(1 — qx/p)(\ — rx/p) ' 



0(0,,)- QM 



'1 -py/q)(l - ry/q) 



) \/l — xqrw 2 — — (1 ) \/l — prw' 

pw ) p \ qw J 

( ( 1 — ) a/1 — yprw 2 — — f 1 ) \/l — qrw 2 ) 

\\ QwJ q V pwj J 



where w is the smallest positive solution ofw = 2+pqrw 3 . The complete generating function Q(x, y) 
can be obtained using (|13|) . When p = q, then w = 1/p, and the above expressions simplify to 

Q(x,0)=Q(0,x)- Q( °' 0) 



■y/l — rx/p 

If r > mm(p,q), no solution of converges on < 1. 

Observe that exchanging p and q, and x and y, leaves the solution unchanged, in conformity with 
the diagonal symmetry of the model. From the algebraic equation defining w, it is not difficult to 
see that 

(pw — 1) 2 (1 — qrw 2 ) = (qw — 1) 2 (1 — prw 2 ) = (rw — 1) 2 (1 — pqw 2 ). 
Moreover, an elementary study of the function f(z) = pqrz 3 — z + 2 gives bounds for w: 

1111 

T7 < W < —— < — < -, (15) 
M Jpq m r 



where m = min(p, q) and M = max(p, q), with equalities holding if and only if p = q. Hence 



< (Mw — 1) yl — mrw 2 = —(mw — 1) yl — Mrw 2 = —(rw — — pqw 2 , (16) 

and this allows us to rewrite the second part of the expressions of Q(x,0) and Q(0,x) in various 
ways. This will be useful in Section 13.31 where we make further comments on this solution, and 
relate it to Flatto and Hahn's formulation. We shall prove Theorem 0] in Section 13.21 For the 
moment, let us derive from it the stationary distribution of the Markov chain of Figure |21 This 
result is actually not given explicitly in |14j . 

Corollary 5 (The stationary distribution) Assume r < min(p, q). Let w be the smallest pos- 
itive solution of w = 2 + pqrw 3 . The Markov chain schematized in Figure |2J with the additional 
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condition ((12)). has a unique stationary distribution (pij), given by 

w(v — r)(q — r)\p — q\ 

Po,o 



6pq(l — rw)^Jl — pqw 2 
^Tpifix' = ^(Q(x,0)-Q(0,0)), 



i>0 



j>o r 

P^ xi y j = -(Q(x,y)-Q(x,o)-Q(o, y ) + Q(o,o)) 



i>0,j>0 

where Q(x,y) is the function of Theorem® taken with Q(0,0) = poo- When p = q, then the 
expression o/j>o,o should be taken to be 

P0fl = \{l-r/pf' 2 . 
If r > mm(p, q), then the Markov chain has no stationary distribution. 

Proof. The stationary distribution is related to the series Q{x,y) satisfying 1)13)1 by 1)14)1 . so that 
the expressions of the three series above are obvious. We only have to determine which value of 
Q(0,0) guarantees the normalizing condition ()!!)) . This condition reads 

Q(0, 0) + ~(Q(1, 0) - Q(0, 0)) + 1(Q(0, 1) - Q(0, 0)) + -(Q(l, 1) - Q(l, 0) - Q(0, 1) + Q(0, 0)) = 1, 
that is, 

Q(l,l)-gQ(l,0)-pQ(0,l)=r. (17) 
When y = 1, a factor (1 — x) comes out of ((TB*)) . leaving 

(p-rx)Q(s,l)=pQ(0,l). (18) 

Setting x = 1 gives 

Q(l,l) = — Q(0,1), 
p — r 

and of course, a symmetric argument yields 

Q(l,l) = -^-Q(l,0). 
g — r 



Hence the normalizing condition 1)17)1 reads 



Q(1,0) = £-^. (19) 
3<? 



Now, from the expression of Q(x,0) given in Theorem 0J we obtain, if p ^ q, 

Q(l, 0) = — pQi^fy ((pui — 1) y/l — qrw 2 — (qw — 1) y/l — prw 2 \ . 

wyp — q)(p — r) V / 
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Using H16|) . the above expression for (5(1,0) can be rewritten as 



_ 2pQ(0,0)(l-r W ) v /l-pg W 2 
— — rj 

and the condition l|19j) gives the value of Q(0,0) = po,o- 

When p = q, the simplified expression of Q(x, 0), given in Theorem 01 gives Q(1,0) 
Q(0, 0)/yl — r/p, and the result follows. 



3.2 Proof of Theorem H 

Our solution of Eq. (|13|) follows the same idea as Section r2.3l we shall exploit an invariance property 
of the kernel. However, we do not have the length variable t any more, which means that we are 
no longer in a power series context, but rather in the world of functions of two complex variables 
x and y. Consequently, certain operations that were performed formally in Section 12.31 fe.g.. the 
extraction of coefficients) now need to be justified analytically. The analytic lemmas we need are 
gathered in Section T3.2. II below. In Section T3.2.21 our main functional equation (|13|) is transformed 
into two functional equations defining two functions U(x,y) and F(x,y), which are respectively 
symmetric and anti-symmetric in x and y. Section [3.2.31 is the heart of the proof: there we solve 
these two equations, using the algebraic kernel method. The reader may skip directly to the latter 
section in order to recognize the logic of the kernel method, transferred to an analytic context. 



3.2.1 Preliminary results 

Our first lemma tells us that the domain of convergence of Q(x, y) is actually larger than the unit 
polydisc |x| < 1, \y\ < 1. 

Lemma 6 Let Q(x,y) be a power series solution of (|13|) that has nonnegative coefficients and 
converges for \x\ < 1 and \y\ < 1. Then this series also converges in the following domain: 

{\x\ < p/r, \y\ < 1} U {|x| < 1, \y\ < q/r}. 

Moreover, 

Q(0,1) Q(1,0) 
Q(x, 1) = t- and Q{1, y) - 



1 — rx/p 1 — ry/q 

In particular, such a solution Q(x,y) can only exist if r < min(p, q): if this inequality does not 
hold, the Markov chain has no stationary distribution. 

Proof. For \x\ < 1, the expression of Q(x,l) follows from (|18|). The expression of Q(l,y) is of 
course symmetric. These two series must converge when x = 1 and y = 1: this forces r to be 
smaller than p and q. Given the nonnegativity of the coefficients, the values of Q(x, 1) and Q(l, y) 
imply the convergence of Q(x,y) in the desired domain. 

■ 

The extraction of coefficients will be based on the following result |31l Chap. 10, Ex. 25]. 
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Proposition 7 Let f(z) be an analytic function in the annulus A = {r < \z\ < R}. There exists a 
unique bi-infinite sequence (a n ) ne % such that for all z € A, 

f(z) = J2anZ n . 

neZ 

Moreover, the convergence is absolute. In other words, f is the sum of a function analytic in the 
disk {\z\ < R} and a function analytic for \z\ > r. 

It will be convenient to work with an equation whose kernel is symmetric in x and y. In Eq. (|13|) . 
let us replace x by px and y by qy. Multiplying through by xy gives 

(xy - x - y - pqrx 2 y 2 )Q(px, qy) + x(l - qy)Q(px, 0) + y(l - px)Q(0, qy) = 0. (20) 

Let K(x, y) = xy — x — y — pqrx 2 y 2 denote the kernel of this equation. The discriminant of K, 
taken as a polynomial in y, is (x — l) 2 — Apqrx s . Let us denote 

A(x) = (1 — x) 2 — Apqrx. 



Lemma 8 Assume r < min(p, q). Let m = min(p, q) and M = max(p, q). The three roots of A(x), 
denoted Xi, i = 0,1,2, are real, and satisfy 

1 1 1 

< Xq < 1 < x\ < — < — < - < x 2 - 
M m r 

Proof. The variations of A(x) are easy to study. Note that A(x) is a square when x = 1/p, 1/q, 
or 1/r. 

■ 

Recall that w is defined as the smallest positive solution of w = 2 -\-pqrw 3 . This implies that 
l/(pqrw 2 ) = X2- The lemma above gives the following factorization of A, which will play the role 
of the canonical factorization of Section [2 

A(x) = A A + (a;)A_(x) (21) 

with 

Ao = 4pqrx2 = 4:/w 2 , A + (x) = 1 — x/x2 = 1 —pqrw 2 x, A_(x) = (1 — xxq)(1 — xx±). (22) 
As a polynomial in y, the kernel K(x, y) of Eq. ()20|) has two roots: 

Zpqrx Zpqrx 
The elementary symmetric functions of the are polynomials in x = 1/x: 

Y + Y 1 = ~ 2) and Y Y 1 = —. 

pqr pqr 

Using the canonical factorization of A, we see that Y$ and Y\ are at least analytic in the annulus 
x\ < \x\ < X2- Let us study these functions a bit more precisely when x is real. The following 
lemma is illustrated by Figure ^] 
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Figure 4: The real branches of the functions Yi, for p = 1/3, q = 1/2 and r = 1/6. 



Lemma 9 We still assume that r < mm(p, q). The functions Yq and Y\ are well-defined and real 
for x E (— oo, xq] U [x\, 22]. In particular, 



Y {l/M) = l/m, Y (l/m) = Y (l/r) = 1/M, Y x (l/M) = *i(l/m) = 1/r, Yi(l/r) = 1/m. 



Each of the derivatives Yq(x) and Y{(x) admits a unique zero on the interval [xi,X2], respectively 
denoted by V2 and v\ . Moreover, 



The function Yq decreases between x\ and V2, and increases between V2 and X2, while the function 
Y± increases between x\ and v\, and then decreases up to X2- 
Finally, for x G (l/m, ^2), one has 



Proof. The proof is a bit tedious, but elementary. We merely sketch the different steps. 

The first assertion comes from the study of the discriminant A(x) (Lemma IB}- The values of 
Yi at the points 1/M, l/m and 1/r are obtained by a direct calculation. 

Let us now focus on the interval [xi,X2]- Given that (1 — x)Yi = 1 + pqrxY? and x > x\ > 1, 
we have < Yq(x) < Y\{x). The derivatives of Yq and Y\ with respect to x can be written: 



Given that Yq and Y\ are positive on [£1,22], anv ro °t °f these derivatives will be larger than 2. 
The equation satisfied by the Yi implies that these roots are also solutions of x — 2 = pqrx 3 . The 
polynomial pqrz 3 — z + 2 has two roots larger than 2. Let us denote them v\ and t>2, with v\ < V2- 
Note that v\ is actually the number w defined in Theorem QJ If x = Vi, then A(x) = (x — 3) 2 . 
Hence Vi belongs to the interval [21,22]. Evaluating the numerators of Yq and Y{ at v\ and V2 



X\ < — < Vl = W < — < V2 < - < X2- 

M m r 
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shows that v\ cancels Y{, while V2 cancels Yq. Finally, we compute the numerators of Y ' and Y{ at 
x\ and X2- This determines the sign of these derivatives and completes the study of the variations 
of Yq and Y\ . 

The last assertion is proved by studying the function y t— > K(x,y), for x fixed. 

■ 

Let x € (xi,x%) and let y belong to the annulus {Y (x) < \y\ < Yi(x)}. Let K r = K/(xy) = 
1 — x — y — pqrxy be the rational version of the kernel. Then the following expansion is convergent: 

k =vk> G=k + - =7^ + g ™j • (23) 

Let F(y) be analytic in the same annulus. Let us write 

F(y) = M = F ~(y) + /o + F+ (y)> 

where .F + and F~ are the positive and negative parts of F. Then 



= _J^( F -(i /yi) + /o + j p+(y o) ) 
A r y/A{x) 

/ F-(pqrxY ) + f +F+(Y )) 



(24) 



In particular, if F(y) = G(y) — G(xy/(pqr)), then 

[y°]^ = o. (25) 

3.2.2 Simplification of the functional equation 

Let us go back to Eq. (f2*0|) . The function Q(x, y) = Qipx, qy)/((l—px)(l — qy)) satisfies an equation 
that is symmetric in x and y: 

(xy - x -y - pqrx 2 y 2 )Q(x, y) + xQ(x, 0) + yQ(0, y) = 0. (26) 

Yet, we shall see that this function is, in general, not symmetric in x and y (Section I3.3|) . We 
shall symmetrize it by considering Q(x,y) + Q(y,x). More precisely, we shall study separately the 
functions S(x,y) and D(x,y) (as in Sum and Difference) defined by 

S(x,y) = (l-qx)(l-py)Q(px,qy) + (l-px)(l-qy)Q(py,qx), 
D(x,y) = {l-qx)(l-py)Q(px,qy)-(l-px)(l-qy)Q(py,qx). 

By Lemma these functions are analytic in 

V={\x\< 1 -,\y\<^}u{\x\<±\y\<l} (28) 
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with M = m&x(p,q). Note that S(x,y) and D(x,y) satisfy the same equation: 

S(x, y)(xy — x — y — pqrx 2 y 2 ) + x(l — py)(l — qy)S(x, 0) + y{l — px){l — qx)S(0, y) = 0, 
D(x, y)(xy — x — y — pqrx 2 y 2 ) + x(l — py)(l — qy)D(x, 0) + y(l — px){l — qx)D(0, y) = 0. 

(29) 

However, S(x,y) is symmetric in x and y, while D(x,y) = —D(y,x). We shall solve separately the 
two equations, taking into account the respective symmetry or anti-symmetry condition. For each 
equation, we will obtain a unique solution, up to a multiplicative factor. A relation between the 
two factors will be found by noticing that, by definition of S and D, 

S( -,o] +D[ -,o] =0. (30) 



.9 

First we rewrite the above equations on S and D, by expressing x(l — py)(l — qy) in terms of 
the kernel: indeed, 

rx(l — py)(l — qy) = (r — x)(x + y — xy(p + q)) — xK(x, y). 
Let K r = K/(xy) = 1 — x — y — pqrxy. The equation satisfied by S can be rewritten as 
xyrS(x,y) -yS(x,0) -xS(0,y) 

x + y-xy(p + q) (dl) 

= xT(x) + yT(y) 

T{x) = (1 - rx)S(x, 0) = (1 - rx)S(0, x). (32) 



with 
Similarly, 



xyrD(x,y)-yD(x,0) -xD(0,y) nU ,- , n , n v 

K r ; t — ■ — > = (x - r)D(x,0) + (y - r)D(0,y) 

x + y — xy(p + qj 

= E(x)-E(y), 



where 



E(x) = {x- r)D(x, 0) = -(x - r)D(0, x). (33) 

We have taken into account the fact that D(x, y) is anti-symmetric, so that, in particular, D(0, 0) = 
0. The functions T and E are analytic for \x\ < 1/r. 

From the fact that the curve x + y = xy(p + q) intersects the domain of convergence of S(x, y) 
(near the origin), we derive from 1)31 J) the existence of a function U(x,y), analytic in T>, such that 

xyrS(x, y) - yS(x, 0) - xS(0, y) = (x + y - xy(p + q))U(x, y). 

A similar statement holds for the function xyrD(x,y) — yD(x,0) — xD(0,y). But this function 
vanishes as soon as x = or y = 0, so that we can actually write 

xyrD(x, y) - yD(x, 0) - xD(0, y) = xy(x + y - xy(p + q))F(x, y), 

for a function F that is analytic is T>. Finally, we shall need the following initial conditions 

U(x, 0) = -5(0, 0) = -2Q(0, 0), T(0) = 2Q(0, 0). (34) 

The equations ()29|) have thus been replaced by the following simpler equations: 

K r U(x,y) = xT(x) + yT(y), (35) 

K r xyF{x,y) = E(x) - E(y), (36) 

where K r = 1 — x — y — pqrxy and U, T, F, E are analytic in the domain D defined by (|28|) . 
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3.2.3 The algebraic kernel method 

We apply the algebraic kernel method of Section 12.31 to the equations (jSMSSl) • We observe that K r 
satisfies the following invariance condition: 

K r (x, y) = K r ( — , y ) = K r ( x, — ) = K r . 
\pqr J \ pqrj 

Recall that m = min(p, q) and M = max(p, q). Let us fix x in the interval (1/m, 1/r) and restrict 
y to the annulus 

Yo(x) < \y\ < ±. (37) 
By Lemma El this annulus is non-empty, and moreover, 

— < \y\ < Y x {x). (38) 
pqx 

The pairs (x,y) and (xy/(pqr),y) both belong to the domain of convergence T>, and we thus have, 
in addition to (|35h36|) : 

K r U[^-,y) = pqrxy T ) + yT(y), 

\pqr J \PQf , 

2 - ^M. ,\ r ( xy\ 



K r ——F ( —^-,y ) = -E(y). (40) 

pqr \pqr J \pqr J 



A linear combination of l|35j) and (|39|) gives 



2U(x,y) 



( xy \ 1 / / xy 

U ,y I = I 2xT(x) + yT(y) - pqrxyT 

\pqr J K r \ \pqr 



In view of (|37H38|) . the expansion of 1/K r given by (|23|) is convergent. Recall that x is fixed; we 
can now use (|25[) to extract from the above equation the coefficient of y°. We obtain 

^(..OJ-^f-i-)-^ (41) 
\pqrj y/A(x) 

where U^{x) denotes the diagonal of the series U(x,y): 

U(x,y) = ^2 u itj x l y j =>• U d (x) = ^u^x 1 . 

i,j>0 i>0 

Given that U(x,y) converges absolutely in the domain T> given by (|28|). the sub-series Ud(x) is 
convergent for \x\ < l/(rM). By analytic continuation, (|41|) holds in the annulus {1/m < \x\ < 
1/r}. Recall that U(x,0) is actually a constant (see (|34(l ). We now use the canonical factorization 
of A(x), given by (f2l|). We multiply (jUJ by y / A_(S): 



VAl(i) ( 4Q(0, 0) + ^ f ±) ) = -2=^- 
V \PqrJJ ^/A A + (x) 
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Using Proposition [3 we can extract the nonnegative part of this function. Given that Ud(0) 
£7(0,0) = -2Q(0,0) we obtain, using J22J): 



T(x) = 2Q(0, 0) (1 - TaTM- (42) 

In view of ()31I32|) . we have completed the determination of the Sum function S(x,y). 

Let us now work with the equations (|36|) and (|40|) . Let us divide (|36|) by K r and extract the 
coefficient of y°. We obtain 

E{x) = E(Y Q ). 

If we do the same with (|40[). we simply find = (recall that y) is antisymmetric). However, 
if we extract instead the coefficient of y, we obtain, using (124)): 



^Frf — ) =£(0)-£(Y ), 
pgr ypqr y 

where 

F(x, y ) = ki xi y j => *l0«0 = 5^/^+1^. 

ij>0 i>0 

By combining both equations, and extracting the nonnegative part, one sees that E{x) is actually 
a constant E(0). In view of l)32H33)) and 1)42 jl . one has 

5(x, 0) = 5(0, x) = 2^^ ^ _ ^ V / A^(x) and D(x, 0) = -D(0, x) - 



. /v r v / — \ 7 / — \ ' / 1 

1 — rx \ wJ 1 — rx 

The identity l)3U|) completes the determination of E 1 : 



£(0) = -2Q(0, 0) U - ~) VA+(l/g). 



We can now express Q(x,0) and Q(0, x) explicitly, using (|27j). Thanks to (|16[> . this gives exactly 
Theorem |1J 



3.3 Comments on the solution 
3.3.1 Asymptotics 

For a good understanding of the solution of Theorem 21 or, equivalently, of the stationary distribu- 
tion of Corollary[51 it is useful to determine the dominant singularities of the functions Q(x,0) and 
Q(0, y), and hence the asymptotic behaviour of the numbers ptfi and poj. This is why we briefly 
rederive below a result already proven in [18] . 

Proposition 10 Assume r < min(p, q). The asymptotic decay of the stationary probabilities pi$ 
depends on the relative values of p and q: 

• If p = q, then Q(x,0) is the reciprocal of a square root, and, as i goes to infinity, 

Pi,0 ~ c (r/pf i~ 1/2 

for some positive constant c. 
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• If p < q, then Q(x,0) has a simple pole atp/r as its unique dominant singularity. The decay 
of the numbers pi^ is given by 

Pi,o ~ c {r/p)\ 

• If p > q then Q(x,0) has a square root singularity at l/(qrw 2 ) as its unique dominant singu- 
larity, and 

Pifi ~ c (qrw 2 ) l i~ 3 / 2 . 

Proof. We use standard results that relate the singularities of a series to the asymptotic behaviour 
of its coefficients (see, e.g., [T7j). 

When p = q, the result is clear in view of Theorem^ Otherwise, the three possible singularities 
of Q(x,0) are p/q,p/r and l/(qrw 2 ). The inequalities (fTH|) imply 

p p 1 
- < - < 



r qrw 2 



Hence our first candidate for the radius of Q(x,0) is p/q. However, the numerator of Q(x,0) 
vanishes at this point, so that there is no pole at p/q. Our next candidate is p/r. For this value of 
x, the numerator of Q(x,0) is 

1 A /i — o Q (, 1 



1 ] yl — pqw 2 (1 ] yl — prw 2 . 

rw J r \ qw ) 

According to ()16j) . the first term in this difference is negative. If p < q, then the second term is 
positive. Hence the difference is negative, and Q(x,0) has indeed a simple pole at p/r. 

However, if p > q, then ((TB|l shows that the numerator of Q(x,0) cancels at x = p/r, so that 
the only singularity of Q(x,0) is a square root singularity at l/(qrw 2 ). 

Note that one can compute explicitly, in the same way, the multiplicative constants denoted c 
in the proposition. 



3.3.2 An asymmetry of the solution 

As observed above, the function Q(x,y) = Q(px,qy)/((1 — px)(l — qy)) satisfies an equation that 
is symmetric in x and y (Eq. (J2HJ)). Hence we could expect Q(x,y) to be a symmetric function of 
x and y. This is equivalent to the condition 

Q(x,0)-Q(0,x) = ^)_^1 = . 

1 — px 1 — qx 

However, we derive from Theorem |1] and Eq. (jl6|) that 

Q(x,0)-Q(0,x) 



2Q(0,0)x(pw -l)y/T — qrw* 



w(l — px)(l — qx)(l — rx) 

By (|15[). this quantity differs from 0, unless pw = 1, which forces p = q. Hp = q, the solution satisfies 
Q(x,0) = Q(0,x), so that the symmetry property naturally holds. Otherwise, the asymmetry of 
the result comes from the asymmetric conditions we have required: Q(x,y) must converge when 
\px\ < 1 and \qy\ < 1. 
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3.3.3 Flatto and Hahn's expression 

Assume p < q. Eq. Q16J1 shows that that the numerator of Q(0,y) vanishes when y = q/p and 
V = ( ?/ r - Let us denote 8{y) = yl — yprw 2 . Then the numerator of Q(0, y) is a polynomial in 5(y), 
of degree 3, and two of its roots are 5(q/p) = yl — qrw 2 and S(q/r) = \J\ — pqw 2 . The third root 
is then easily determined to be — yl — prw 2 . Hence, up to a multiplicative constant independent 
of y, the numerator of Q(0,y) factors as 



1 — yprw 2 — yl~ — gru; 2 ) ( yl — yprw 2 — \/Y^pqw 2 ) ( y/l — yprw 2 + \/l — prw 2 



The denominator of Q(0, y) is already factored in y, and also vanishes at y = q/p and y = q/r. Up 
to a multiplicative constant, it factors as 



( y 7 ! — yprw 2 + \/l — qrw 2 ^ ( \/l — yprw 2 — \/l — qrw 2 ^ 



yl — yprw 2 + \J\ — pqw 2 ) ( y 7 ! — yprw 2 — \/l — pqw 2 



Two simplifications occur, and finally 



Q(0,y) = Q(0,0)||| (43) 



where 

T . . \/l — yprw 2 + yl — prw 2 



(-^/l — yprw 2 + y 7 ! — q , ra 2 )(y / l — yprw 2 + yl — pqw 2 ) 
Now, using ()16|) . the function Q(a;,0) can be rewritten as 

Q(x, 0) = — t~T7T" rr if 1 — ~~ ) V 7 ! - xqrw 2 + x ( 1 —J y 7 ! - grw 2 

(1 — qx/p){l — rx/p) \ \ pw J \ pw J 

In this form, the numerator of Q(x, 0) now looks more like the numerator of Q(0, y). More precisely, 
denoting the latter numerator by P(5(y)), the former numerator is exactly — P(— yl — xqrw 2 ), and 
hence factors as 



1 — xqrw 2 + yl — qrw 2 ) ( yl — xqrw 2 + yl — pqw 2 ) ( yl — xqrw 2 — \/l — prw 2 



The denominator of Q(x,0) is also easily factored; two simplifications occur again, and we end up 
with 

Q(x,0)=Q(0,0)||| (44) 

where 



v/l — xqrw 2 + yl — qru; 2 



(yl — xqrw 2 + yl — pra 2 )(y / l — xqrw 2 — yl — pqw 2 ) 

The expressions Q43|) and ()44|1 are, with our notation, the forms given in Flatto and Hahn's pa- 
per [1S| . They are nicely factored, and it is easy to derive from them the singularities of Q(x,0) 
and Q(0,y). However, they have two drawbacks: first, they are only valid when p < q, and hide 
the symmetry of the result in p and q, which is clear from the expressions of Theorem @J Secondly, 
they somehow contain "two many" radicals, and suggest that Q(x,0) and Q(0,y) will be algebraic 
of degree 3 x 2 4 over the field Q(p,q,x,y), whereas as suggested by Theorem HJ they have only 
degree 3 x 2 2 = 12. This can be checked using a computer algebra package, like Maple. 
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4 Enumeration &; probability: the law of the chain 



In this section, we consider again the Markov chain illustrated in Figure 01 We start this chain at 
time at the origin of the lattice, and address the question of computing the probability pij(n) 
that the walk reaches the point at time n. This question is, in essence, close to Section 2: we 
are again enumerating paths according to a certain weight. This weight is the probability that the 
trajectory begins with this path. But this question is also related to Section 3, since we expect the 
probability pij(3n — i—j) to converge to 3pi as n goes to infinity, when r < min(j), q), where pij 
is the stationary distribution of the chain (the factor 3 accounts for the periodicity of the chain) . 

The notation we adopt is similar to that of Section 3: we introduce the following four generating 
functions for the probabilities pij{n): 

Po,o = J>o,oMi n , Pi(x) = M 71 )^*". Mv) = E poj W* n . 

n>0 n,i>0 n,j>0 

n,i,j>0 

The step by step construction of the walks gives the following functional equation: 

(1 — pxt — qyt — rxyt)P(x, y) + (1 —pxt — r' xyt)P\(x) + (1 — q"yt — r"xyt)P2(y) + (1 — xyt)Po t o = 1. 

Again, we assume that the transition probabilities on the border of the quadrant are related to 
those inside the quadrant by the conditions (|12|). This allows us to rewrite the above functional as 

(1 — pxt — qyt — rxyt)Q(x, y) + q(yt — l)Q(x, 0) + p(xt — 1)Q(0, y) = r 

with 

Q(x,y) = P ,o + r'P 1 (x) +r"P 2 (y) +rP(x,y). (45) 
It will be convenient to have a kernel symmetric in x and y, and our starting point will actually be 

(xy — t(x + y + pqrx 2 y 2 ))Q(px, qy) + (i — qy)xQ(px, 0) + (t — px)yQ(0, qy) = rxy. (46) 

We are back to the (safe) world of formal power series in t with coefficients in Q(x, y), and we will 
mimic the obstinate kernel method of Sections 12.11 and 12.21 The only new difficulty arises from the 
absence of symmetry, since Q(x, 0) / Q(0, x) when p ^ q. 

The kernel of the above equation, considered as a polynomial in y, has two roots, 



_ r . . 1 - tx — J(\ — tx) 2 — 4pqrt 2 x 9 , 

Y (x) = - 1 — = t + xt 2 + 0{t- 

Ipqrtx 



1 — tx + y(l — tx) 2 — 4pqrt 2 x 



YAx) = ^ - ' v ^ = _ t-xt 2 + OU s ). 

2pqrtx pqrt pqr 

The elementary symmetric functions of the Yi are again polynomials in 1/x: 

Y + Y 1 = ^—p and r *i = — ■ 
pqrt pqr 



2 
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The discriminant A(x) = (1 — tx) 2 — 4pqrt 2 x vanishes for three values of x: two of them, say Xq 
and Xi, are power series in \ft while the third one, X 2 , is a Laurent series in t that starts with a 
term in t~ 2 . Let us define Z = Z(t) to be the unique power series in t such that 

Z = 1 + 4pqrt 3 Z 3 . 

Then ^pqrt^X^Z 2 = 1, and the canonical factorization of A(x) reads 

A(x) = A A + (x)A_(x) 

with 

A = Apqrt 2 X 2 = —;, A+(x) = 1 - x/X 2 = 1 - 4pqrt 2 Z 2 x, (47) 
Z z 

A_(x) = (1 - xl )(l - xXi) = 1 - tZ(l + Z)x + t 2 Z 2 x 2 . (48) 

As in Section it will be convenient to handle two functions S(x,y) and D(x,y), which are 
respectively symmetric and antisymmetric in x and y. We define them by 

S(x,y) = (t-qx)(t-py)Q(px,qy) + (t-px)(t-qy)Q(py,qx), ^ 
D(x,y) = (t- qx)(t-py)Q(px,qy) - (t-px)(t- qy)Q(py,qx). 

Then 

t (xy -t(x + y + pqrx 2 y 2 )) S(x, y) + 

(t - py){t - qy)xS(x, 0) + (i - px)(t - gx)yS(0, y) = G(x, y) + G(y, x), (50) 
t {xy -t(x + y + pqrx 2 y 2 )) D(x, y)+ 

(t - py)(t - qy)xD{x, 0) + (t - px)(t - gx)yL»(0, y) = y) - G(y, x), (51) 

where 

G(x, y) = rxyt(t — qx)(t — py). 
4.1 Statement of the results 

After all the algebraic series we have met, one might expect the probability generating function of 
the law of the chain to be algebraic again. This is, however, only true if p = q. 

Theorem 11 (The symmetric case) Assume p = q. The three-variate generating function for 
the probabilities pij(n) is algebraic, and can be expressed explicitly in terms of the unique power 
series Z = Z{t) satisfying Z = l + 4pqrt 3 Z 3 . In particular, the generating function of walks ending 
at the origin is algebraic of degree 6: 

* - E,,„(3^ . L (VES _ {& _ ,) , 

where A_(x) is given by (|48|) ■ More generally, the series Q(px,0) = Pq,o + r'P\(px) is given by: 



, {2tZ-x)yJ&-{p/t)yJ& + {x) \ 

(t - x{l -p)+ prx t )Q{px, 0) = — I £ _ 2 ^ 2t + x{l - p) I 

where A + (x) is given by (|47|) . The expression of Po,o can be recovered from the value of Q(px,0) 
by setting x = 0. 
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This theorem, and all the results of this section, will be proved in Section 14.21 

What happens in the general case? We have expressed the series Q(px, qy) in terms of two series 
S(x,y) and D(x,y), which are respectively symmetric and antisymmetric in x and y. It turns out 
that the Sum series S(x, y) is always algebraic, while the Difference series D(x, y) is transcendental 
(unless p = q). The algebraicity of S(x, y) has an interesting consequence: The generating function 
-Po.o that counts walks ending at the origin is always algebraic, even when p ^ q. 

Theorem 12 (The general case: algebraic part) The series S(x,y) defined by (|49|) is alge- 
braic and can be expressed explicitly in terms of the unique power series Z = Z(t) satisfying 
Z = 1 + Apqrt 3 Z 3 . In particular, the coefficient of x°y° in S(x,y) is an algebraic series in t. It is 
equal to 2i 2 Po,o> where Po,o counts walks ending at the origin, and we have: 

2pg-P ,o + r(l - r) = A m + A q>p 

where A pq is the following algebraic series in t: 



(p(l-2p)-qrt 3 )^A-(p/t) 
p ' q (l-t 3 )(l-2pZ) 

and A-(x) is given by (|48j) . The algebraic series Po,o has degree 6 if p = q, and degree 12 otherwise. 
More generally, the series S(x, y) satisfies: 

(t-(l- p)x + t 2 qrx 2 ) (t - (1 - q)x + t 2 prx 2 ) ^l + r _B^l = 
v ' ' t 2pq 



' Ap qFp q(x) + AqpFqp(x) 



2pqZ 

where A+(x) is given by (|47|). and F Pt q(x) and H{x) denote the following polynomials in t and x: 

Fp,q( x ) = {t ~ X( l) (t — (1 — q)x + t 2 prx 2 ) , 
H{x) = q(2t-x+px)Fp tq (x)+p(2t-x+qx)F qtP (x)-(p-q) 2 x 2 (2rt-(l-p)(l-q)x+2pqrx 2 t 2 ). (52) 



The expression of Po,o can be recovered from the value of S(x, 0) by setting x = 0. An expression 
for S(x,y) can be obtained using (|5Uj). 

This theorem will allow us to complete the proof of the following result, announced in Section [31 

Corollary 13 The Markov chain schematized in Figure with the border conditions ()12|) . is 
ergodic (that is, has a stationary distribution) if and only if r < min(p, q). 

Theorem^] specializes to Theorem 1111 when p = q. It states that walks ending at the origin have 
an algebraic generating function. What about the generating functions P\{x) and P2(y) that count 
walks ending on the x- or y-axis? By symmetry of the model, P\{x) is algebraic if and only if P^y) 
is algebraic too. In view of Eqs. (|45|). (|46j) and (|49|). this holds if and only if S(x, 0) and D(x, 0) are 
algebraic. If p = q, then D(x,y) is obviously zero, and Theorem II II tells us that all the generating 
functions under consideration are algebraic. If p ^ q, we shall prove that D(x,0) is transcendental 
(but D-finite), and give an explicit expresion of it. 
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So far, we have expressed many of our series in terms of the canonical factorization of the 
discriminant A(x). This is the case, for instance, in Theorem 1121 above, where the expression of 
S(x, 0) involves sj A + (x), which we could call the positive multiplicative part of y/ A(x). In order 
to express D(x,0), we need to introduce the positive additive part of \J A(x), as defined by (J2J). 
More precisely, the expression of D(x,0) will involve the positive (additive) part of 



/A(x) 

B(x) := (Y - Y 1 )(2t - x + pqrtx 3 ) = ^ — — (l - 2tx - pqrtx 2 ) . (53) 

pqrt 

We shall compute below the expansion of B in t and x, using the Lagrange inversion formula. In 
particular, we will see that the positive (additive) part of B{x) reads 

B + (x) = -xt- x 2 + 2C + (x) (54) 

where all terms in C + {x) are multiples of x s . We then define the series C~(x) by 

i _ o^t 

B{x) = + 2C-{x)-xt-x 2 + 2C + {x). (55) 

pqt 

Observe that A(t/p), and hence B(t/p), is a well-defined Laurent series in t. Clearly, C + (t/p) is 
well-defined too: consequently, by difference, we can define C~(p/t) as a Laurent series in t, even 
though it would be meaningless to replace x by p/t in the expansion of C~(x). 

Theorem 14 (The general case: transcendental part) When p ^ q, the series D(x,y) de- 
fined by (|49j) is D-finite but transcendental. The same holds for its specialization D(x,0). Conse- 
quently, the series Pi(x) and P2(y) which count walks ending on the x- or y-axis are transcendental. 
The series D(x,0) satisfies: 

D(x 0) 

it - (1 - p)x + t 2 qrx 2 ) (t - (1 - q)x + t 2 prx 2 ) y ' 1 +x(p-q) (t 2 (l - r)rx 2 - x/2 + 1) = 

rx{p-q)t 2 C + (x) - -J—(pC-(p/t)F M {x) - qC~ {q/t)F q ^{x)) , 
where C + (x) and C~(x) are defined by ()54H55|) and, as in Theorem \\2V 

Fp,q(x) = (t — xq) (t — (1 — q)x + t 2 prx 2 ) . 
An expression of D(x,y) can then be obtained using ()51|) . 

Note the similarities between the expressions of S(x,0) (Theorem 112(1 and D(x,0) (Theorem 114(1 . 
One could take the sum and difference of these expressions to recover the series Q(px,0) and 
Q(0,qx), but, as no significant simplification arises, we shall not do this. 

There is still one natural question that is not answered by the combination of the above two 
theorems: we have seen that the generating function Po,o of walks ending at the origin is algebraic, 
but that the series Pi(x) that counts walks ending on the j;-axis is transcendental. Yet, for i > 0, 
the coefficient of x l in Pi(x), being 

n>0 

counts walks ending at (i,0) and might be algebraic. The following corollary tells us that this is 
not (systematically) the case. 



27 



Corollary 15 Some of the series P^o are transcendental. 

Note. We can obtain an explicit expression of the series C + {x) and C~(x) by expanding B{x) in 
x and t. Let us write 

Y -Y 1 = 2Y - (Yo + Y,) = 2Y - ^ ~ f\ 

pqrt 

and observe that the series Yq, which cancels the kernel, is Lagrangian in t: 

Y = t (l + xY + pqrxY$) . 
The Lagrange inversion formula yields 

V ' pqrt ^ _ k\(k + l)!(n - 2k)\ 



Consequently, 

\n/2\ 



Jk-n+2,sk (3fc - n + l)(3fc - n){n - 2)! 



C+(x) = ^ E fcI ( fc + 1 )!(„ -2*)! 

n>2 fc=[(n-l)/3l V ; V ; 

One may also write an explicit expansion of C~(x). 
4.2 Proofs 

Proof of Theorem 1121 Let us start from the equation (|5U|) defining S(x,y). As in Section [2,11 
the pairs (x,Yq) and (Yq,Yi) cancel the kernel and can be substituted for (x,y) in this equation. 
We thus obtain two equations: 

(t - P Y )(t - qY )T(x) + (t - px)(t - qx)T(Y ) = G(x,Y ) + G(Y ,x), . . 

( t -pY 1 )(t-qY 1 )T(Y ) + (t-pY )(t-qY )T(Y 1 ) = G(Y ,Yt) + G(Y U Y ), [00) 

with T(x) = xS(x, 0). Let us form a symmetric function of Yq and Y\ based on a divided difference: 
We multiply the first equation by 2(t — pY\){t — qY±) and the second one by (t — px)(t — qx), and 
take the difference of the resulting equations: 

2(t_py )(t_gy )(t_ p y 1 )(t- g y 1 )T(x)+(t-px)(t-^)((t-py 1 )(t-( ? y 1 )r(y )-(t-py )(t- (? y )T(y ] 



= 2{t - P Y 1 )(t -qY 1 )[ G(x, Y ) + G(Y Q ,x) )-(t- px)(t - qx) ( G(Y , Y 1 ) + G(Y 1 ,Y , 
Recall that 



K(x, y) = xy — t(x + y + pqrx y ) = —pqrtx (y — Yo)(y — Y\). 

We use this identity to express the coefficient of T{x) as a Laurent polynomial in x and t. Then, 
we separate the symmetric and anti-symmetric parts of the right-hand side using 

*Cro,li) = i(*(y 0) Fi) + ^(Fi,y )) +i(*(Fo,Fi)-*(Fi,y ) 
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This gives 

(t - (1 - p)x + t 2 qrx 2 )(t - (1 - q)x + t 2 prx 2 ) t H(x) 



T'(x)+- „ 

pqr z x^ p z q z rx^ 



(t - px) (t - qx ) I ( t _ pYl )(t - qYi)T(Y ) - (t - pY )(t - gYo)T(Yi) ) + (Y - Y^-^- (57) 



pqx 

where H{x) is given by (|52j) and J{x) is also a polynomial in x and t: 

J(x) = qF p , q (x) +pFg :P (x) + x(p - q) 2 (t - rx). 

As we are getting used to the method, let us merge the next two steps: instead of first dividing by 
(Yq — Yi) and then multiplying by \J A_(x), let us divide (fSTj) by 2 \J A + (x) / (pqrx) = —2t(Yo — 
Y 1 )/ y /A A_(x). We obtain, in view of (07J): 

(t - (1 - p)x + t 2 qrx 2 )(t - (1 - q)x + t 2 prx 2 ) . , t H(x) 

-1 \X)+- 



rx 3 \J A + (x) 2pqx 2 \J A + (x) 

V'^^M f(t- px )(t - qx) {t ~ pYl){t ~ qYl ^ T ^ ~ - gYo)r(Fi) i J(x) 



2Z v t(y -ii) p?z 

The left-hand side of this equation, as a Laurent series in x, has valuation —2, while the right-hand 
side only involves powers of x smaller than or equal to 2. Extracting the positive part in x, and 
mutiplying by x 2 gives 

(t - (1 - p)x + tVz 2 )(t - (1 - g)x + t>x 2 ) + tH(x) = 
rx v / A + (x) 2pq^A + (x) 

where L(x) is a polynomial in x, of degree 4, with coefficients in Q[t, Ti, T2, where Tj = Tj(t) 
denotes the coefficient of x l in T(x). We do not give the explicit expression of L(x), but refer the 
reader to his/her favourite computer algebra system. 

We now have to determine the three unknown functions T\,Ti and T3. Fortunately, we can 
compute T{x) at three values of x using (|56|) : first at x = t/p, then at x = t/q, and finally at 
x = W, where W is the unique power series in t that satisfies K(W, W) = (so that W = Yq(W)). 
Remarkably, W is simply related to the parameter Z by W = 2tZ. The three values of T{x) that 
we obtain are: 

t 3 (p - q) (q-r + ^(1 - p) 2 - 4t 3 qr) 
T ^ ~ 2p»,(l-f) (59> 



t 3 (q — p) (p — r + y(T— q) 2 — 4t 3 pr 



r( ' / "» = wd-f) ' <60> 

T(W) = rtW 



2 



Setting a; = in (|58|) . we find that the left-hand side vanishes. Hence L(W) = 0, and this gives 
an expression of T2 in terms of T\ : 

2r-l 
T 2 = -rt+— Tl . 
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The polynomial L(x) now takes the following form: 



. (W -x)(t(t-x-tpqrx 3 )M + M 1 x 2 ) 

L{x) = p^w (61) 

where 

M =pqT 1 + r(p + q)t 2 , 

and Mi involves both T\ and T$. It remains to evaluate Eq. (|58(1 at x = t/p and x = t/q, using the 
expressions of T(t/p) and T(t/q) given by (|59|) and (|bl)|). to obtain 

_ t 3 / (p(l - 2p) - t*qr)JE-{p/t) (g(l - 2g) - tV)y / A_(gA) \ 
° 1 - t 3 y t-pW t-qW J ' 

t 3 ( {p(l - 2p) - t 3 qr)(q(l - q) + prt 3 ) v r Kjj)/Tj 



Mi 



1 - t 3 V * - pW 

(g(l - 2g) - iV)(p(l -p) + gri 3 )VA_(g/t)' 



Theorem IT2*1 follows, using (|oT]) and (|58jh 



Proof of Corollary 1131 We have already seen that the condition r < min(p, q) is necessary for the 
chain to have a stationary distribution (Lemma EJ). Assume this condition holds. As the chain is 
irreducible, it suffices to prove that the point (0, 0) is positive recurrent, that is, that the probability 
Po t o(3n) converges to a positive constant as n goes to infinity ^Oj. The generating function of these 
numbers, denoted Pofl, is given explicitly in Theorem 1121 

This leads us to determine the smallest singularity of the series A Ptq . The technique is standard 
for algebraic functions, and we only sketch the main steps. The series Z becomes singular at 
t 3 = l/(27pqr) > 1. Then, we note that 

P (i - 2 P) - w e = a-^)(i + (2 P 1)2(1+2^ 

4pZ 6 

so that there is actually no pole in A P)Q if Z reaches l/(2p). Moreover, if A_(p/t) = 0, then 
A.(t/p) = 0. But, as t increases from to 1, A(t/p) = (1 — p) 2 — Aqrt 3 decreases from (1 — p) 2 to 
(q — r) 2 and thus does not vanish. Hence A Ptq (and A q:P ) has its smallest singularity at t 3 = 1, and 
this singularity is a simple pole. Consequently, the coefficient of t 3n in A p>q tends to a constant as 
n — > oo. The same holds for A q ^ p . 

It remains to show that the sum of these two constants is not zero. We are actually going to 
compute them explicitly: this will not only conclude the proof of the corollary, but also allow us 
to recover the value of po,o given in Corollary \5\ First, we note that Z(l) = 1 + ApqrZ(l) 3 , and 
conclude that Z(l) = w/2, where w is the real number defined in Theorem 0J Then, the definition 
of the canonical factorization gives, when t = 1, 

A(l/p) = {q- r) 2 = AoA+(l/p)A_(p), 
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and so by (|I7|). 



2y/l — qrw 2 
Thus, as t — > 1, 

2pqP o ~ W [ (Pi 1 ~ 2 P) ~ g r )(g ~ r ) + ~ 2g) -pr)(p- r) \ 
2(1 — t 3 ) I (1 — pw) \Jl — qrw 2 (1 — qw)\Jl — prw 2 J 

Note that p(l — 2p) — qr = (q — p)(p — r). Using (|16j) . we rewrite the above identity as 

1 w{p — r)(q — r)\p — q\ 



1-t 3 (l-rw)y/T 



pqw 



2 



It follows that, as n — > oo, 

w{p - r)(q - r)\p 

Po,o(3n) 



2pq{\ — rw)^Jl — pqw 2 
Given that the chain has period 3, this agrees with Corollary [SJ 



Proof of Theorem 1141 Let us start from the equation (|5T|) defining D(x,y). The pairs (x,Yq) 
and (Yoi^i) cancel the kernel and can be substituted for (x,y) in this equation. We thus obtain: 

(t-pY )(t-qY )E(x)-(t-px)(t-qx)E(Y ) = G(x,Y ) - G(Y ,x), 
(t-pY 1 )(t-qY 1 )E(Y )-(t-pY )(t-qY )E(Y 1 ) = G(Y , Y t ) - G(Y U Y ), 

with E(x) = xD(x,0) = —xD(0,x). We now want to form a symmetric function of Yq and Y± 
based on a sum. We multiply the first equation by 2(t — pY\){t — qYi) and the second one by 
(t — px)(t — qx), and take the sum of the resulting equations: 

2(t-pY ){t-qY )(t- P Y 1 )(t-qY 1 )E(x)-(t- P x){t-qx)^ 

= 2{t - pY{){t - gYi) (G(x, Yq) - G(Y , x)) + (t - px)(t - qx) (g(Y , Y x ) - G(Y U Y: 

As above, we use the expression of the kernel to express the coefficient of E(x) as a Laurent 
polynomial in x and t, and split the right-hand side into a symmetric and an anti-symmetric part. 
After multiplying by x, we obtain: 

(t - (1 -p)x + t 2 qrx 2 )(t - (1 - q)x + t 2 prx 2 ) t(p - q)I(x) _ 

„„^.2 3 E\X)-\ 2 — 



pqr*x° p^q^rx 



x(t-px)(t - qx)((t-pY 1 )(t - qYi)E(Y ) + (t-pY )(t - qY )E(Y 1 ) ) + (p - q )tl?M. ( 62 ) 



pq 



where B{x) is given by ()53() and 

I(x) = —x 3 pq — t 3 pqrx 3 + pqr(l — 2r)t 2 x 4 + (r + 3pq)x 2 t — (1 + r)xt 2 + t 3 . 
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As a Laurent series in x, the left-hand side of the above identity has valuation —2. The term 
involving E(Yq) and E(Y\) only involves powers of x smaller than or equal to 2. But B(x) is a 
series in t with coefficients in Q[x,x], containing arbitrarily large positive and negative of x, and 
this is where the transcendence of the solution stems from. 

Note that the coefficients of x and x 2 in B(x) are especially simple, being respectively —t and 
— 1. Let Ei = Ei(t) denote the coefficient of x % in E(x). Given that E{x) = xD(x,0) and D(x,y) is 
antisymmetric, we have Eq = E\ = 0. Let us first extract from (j62() the coefficient of x: we obtain 
a relation between E2 and E3: 

T 

E3 = -E 2 + r(q-p). 
Now, extracting the positive part of (|H2*|) and multiplying by x gives 

(t — (1 — v)x + t 2 qrx 2 )(t — (1 — q)x + t 2 prx 2 ) . , , ,B + (i) r , . 
2^ i ^ 2 ^ £ = x (p _ g ) t 3 + ( 63 ) 

pqr z x z pq 

where B + (x) is the positive part of B{x) and is a polynomial in x (of degree 3) with coefficients 
mQ\p,q,t,E 2 ,E 4 }: 

pqrx 3 t 2 + tx — t 2 — r(l — r)x 2 t 3 + x 2 (r 2 — pq) x 2 t 2 

L(x) = -2 7, E2+2 ^ E4 

pqr z pqr z 

,(l-2r)t 3 a . . 2 + 2r + t 3 r 2 . , t 2 

- (p - q)- —x 3 + (p- q)t x 2 -2{p-q) x. (64) 

pq pqr pqr 

We have to determine two unknown functions E% and E4. From (|49|1 and the fact that T{x) = 
xS(x, 0) and E(x) = xD(x, 0), we derive that E(t/p) = T(t/p) and E(t/q) = -T(t/q). We evaluate 
Eq. (|6*3*)) at x = t/p and x = t/q, using the expressions of T(t/p) and T(t/q) given by and (|60|) . 
One thus obtains expressions of S 2 and £4 in terms of A(t/p), y/ A(t/q), B + (t/p) and B + (t/q). 
They become much simpler using 



One finds: 



and 



, \/ A(x) , 9s 1 — 2tx „ . . 

B+(x) = — — 1 - 2tx - pqrtx 2 ) 2C~(x). 

v ; pqrt y Pqt 



2 

E2 = ^{qC-(q/t)-pC-(p/t)) 



E * = 1^3 H 1 ~ ^ " Vt")C-(q/t) - p(l - p - qt 3 )C-(p/t)) - r{P g)( 2r + 1 ) . 

The expression of D(x, 0) = E(x)/x follows from these values, using (|fi~H) and (|54|). 

Let us now discuss the algebraic nature of D(x, 0) (equivalently, of E(x)). If E(x) were algebraic, 
the so would be all the coefficients Ei. In particular, the series L{x) given by (|64j) and occurring in 
the right-hand side of (|63|) would be algebraic too. By difference, the positive part of B(x), denoted 
above by B + (x) would be algebraic, and so would be all its coefficients. But the coefficient of x 3 
in B(x) is 

4 Vt 3fc+2 f» ar) k+1 ^ 

A ^ W kKk + mk + 2 )r 
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As k — > oo, the coefficient of t 3k+2 in this series is asymptotic to (27pqr) k /fc 4 , up to a positive 
multiplicative constant. Because of the factor k~ 4 , this cannot be the asymptotic behaviour of the 
coefficients of an algebraic series ^H^, s ° that our initial hypothesis is false: The series D(x,0) is 
not algebraic. However, the general results on D-finite series recalled at the end of Section 1 imply 
that it is D-finite. 

■ 

Proof of Corollary 1151 Assume all the series -P^o are algebraic. By symmetry, all the series Po,j, 
which count walks ending on the y-axis, are algebraic too. In other words, the coefficient of x % in 
Q(x,0) and Q(0,x) is algebraic. In view of (|49|). this holds for S(x, 0) and D(x,0) as well. 

We already know, by Theorem ll21 that S(x, 0) is algebraic. Let us work with D(x, 0) to obtain a 
contradiction. We thus assume that the coefficient of x % in the series E(x) = xD(x, 0) is algebraic. 
By ()63[) . this implies that the coefficient of x l in B + (x) is algebraic too, for all i. The same 
asymptotic argument as above proves that this is wrong. 
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